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Abstract. Anisotropic flow measurements play a crucial role in understanding the 
physics and bulk properties of the system created in heavy ion collisions. In this talk I 
briefly review the most important results obtained so far, recent developments in the 
analysis techniques and the interpretation of the results, and what should we expect 
next, both at RHIC and LHC. I also discuss event anisotropies sensitive to the strong 
parity violation effects. 

Analysis of the event anisotropies (anisotropic flow) in multiparticle production 
in non-central nuclear collisions appeared to be one of the most informative paths 
in understanding the physics and characterizing the properties of the dense and hot 
strongly interacting medium. It has been observed a continuous increase in the value of 
in-plane elliptic flow (t> 2 > 0) [1] from the top AGS energies to RHIC. At RHIC, strong 
elliptic flow [2] comparable in strength to the predictions of ideal hydrodynamics, and 
the hadronization via quark coalescence following from constituent quark number scaling 
of differential flow [31 H], together with the jet quenching, are the key ingredients in the 
picture of sQGP (strongly coupled Quark Gluon Plasma). 

The field is rapidly developing and evolving. From the "discovery phase" of the 
first years of RHIC operations it has been transforming into detailed quantitative 
description of the sQGP phase and subsequent hadronization. The plots of v^je vs 
particle density [31 El E] that have been extensively used in the assessment of the level 
of thermalization reached in the relativistic nuclear collisions and applicability of ideal 
hydrodynamics comes under scrutiny: are the measurements of elliptic flow precise 
enough, are the anisotropies what we think they are and how much are they modified 
by fluctuation processes? When comparing to hydrodynamical calculations, are the 
proper initial conditions used in the calculations? Could it be that we have missed some 
important physics building the models? In the last couple years significant progress has 
been reached answering each and every of these questions. The role of viscosity, flow 
fluctuations, initial conditions (eccentricity and initial flow (velocity) field) are a few 
questions to report in this review. I also discuss preliminary measurements of azimuthal 
"out-of-plane" anisotropies that are related to one of the fundamental questions of the 
strong parity violation. Finally, I briefly overview future measurements at RHIC and 
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LHC. Unfortunately, the space does not allow any detailed discussion of many other 
very important developments, such as flow of 0-mesons [JJ |8] that is important to 
understand the relative flow development in the partonic and hadronic phases, flow of 
the deuterons [7J and if*-mesons as tests of coalescence, heavy flavor flow, etc. The 
KE scaling of elliptic flow observed by PHENIX Collaboration [5] is still not fully 
understood/appreciated. There is no doubt that these measurements will enrich our 
understanding of the dense QCD medium even further. 

Taking into account the significance of the t^/e plot in establishing the sQGP 
picture, I spend most of my time discussing recent developments related to this plot. 
There have been a number of important findings: (a) Along with several indirect 
indications that even in central Au+Au collisions the thermal equilibrium is not 
complete, it was found that even very small, compared to the conjectured low limit of 
r]/s values, the viscous effects lead to a significant reduction in the predicted elliptic flow 
compared to the ideal hydro case. This would lead to a contradiction with experimental 
measurements if other effects, responsible for an increase of elliptic flow (compared to 
the "standard" hydrodynamical calculation) would not be identified. Several of such 
effects have been reported, (b) First, it was noticed that ideal hydro calculations, if 
tuned to describe spectra, yield larger elliptic flow than thought previously, (c) It was 
shown, that in some models, e.g. CGC, the initial eccentricity can take significantly 
larger values than in optical Glauber model that is usually used in hydro calculations. 
The larger eccentricities inevitably lead to larger elliptic flow, (d) Flow fluctuations, 
the nature of which is much better understood in the last years, lead to an increase of 
apparent flow, (e) Finally, it was noticed that the gradients in the initial velocity field 
also lead to the increase in final values of elliptic flow. Any of the above mentioned 
effects can be quite significant, leading each to 20-30% or even larger change in values 
of V2- The final "assembly" of these effects into one reliable model is still under way. 

An attempt of model independent analysis of v%je dependence on particle density 
based on parametrization in terms of Knudsen number has been developed in [TOl [TTJ. 
Using an expression v 2 /s = (v2/£)hydro(l + K/K Q )~ l (where the parameter K rs 0.7 
is independently estimated from comparison to model calculations) to fit the data, see 
Fig. 1, the authors conclude that at RHIC we might be still up to 30% below the ideal 
"hydro limit" even for the most central collisions. Their estimate of the viscosity yields 
values oirj/s = 0.11 — 0.19 depending on the CGC or Glauber initial conditions. Similar 
fits to the STAR data performed by R. Snellings [12] and collaborators lead to similar 
conclusions. 

The magnitude of the viscous effects could be judged already from the early 
calculations [13] where the hydro dynamical evolution at some intermediate stage was 
joined to the transport model to simulate late (viscous) evolution of the system. Later, 
viscosity was attempted to be introduced directly into hydrodynamic calculations [14J. 
Recently there have been performed a few "full" calculations [151 Ell ELY] of the 
hydrodynamical expansion with viscous terms explicitly included into equations. Note 
that the latter (namely the form of these terms) is not that everybody agrees on [16J, 



Anisotropic flow: Achievements, Difficulties, Expectations 



3 



Au-Aii ' — ■ 
Cu-Cu ' — B— 

tit a =4,3 mb 



fir 



0.3 
0.25 
0.2 h 
•J., 0.15 
().[ 
0.05 


0.2 0.4 0.6 0.8 1 1.2 1.4 
(1/SXdN/dyXmb" 1 ] 

Figure 1: Fit to v^/e in terms of Knudsen 
number [TT] . 
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Figure 2: Viscous hydro calculations [15 
compared to data. 



though the difference in the results due to use of somewhat different equations are likely 
small. At the same time everybody agrees on the significance of the viscous effects even 
for the "minimal" values of viscosity (r//s = l/(47r)). The results presented in Figs. [2] 
and [3] show that even the minimal viscosity lead up to ~25-30% reduction in flow values 
in Au+Au collisions and probably more than 50% in Cu+Cu. 

Note that viscosity coefficients calculated in pQCD are usually much larger than 
would be allowed by the data. In this sense, noteworthy are the recent calculations [18J , 
which emphasize the importance of taking into account 2^3 processes. With these 
effects included, the viscosity coefficient appeared to be about an order magnitude 
smaller compared to previous calculations and fall into the "allowable" by the data 
range. 

One can wonder how is it possible that with viscous effects to be that strong that 
the ideal hydrodynamical calculations results are not much higher than the experimental 
measurements? Indeed, there have been identified several effects which possibly lead to 
a significant increase in predicted flow. Taken together with viscous effects they may 
restore the agreement with experiment. 

Firstly, in [19] it was explicitly demonstrated that the ideal hydro calculations can 
not be "tuned" to describe both spectra and elliptic flow. If one tunes the model 
to describe spectra, the elliptic flow values appeared too large, for about 20-30%, 
compared to the data, see Fig. |4| What leaves even more space for viscous effects, 
is the observation [201 EI] that the initial eccentricity calculated in the CGC model 
yields values up to 50% larger compared to the "standard" optical Glauber model. The 
effect has been further studied in [22] ; in [23] it was discussed how eccentricity depends 
on details of CGC model, which can be taken as possibility to investigate CGC model 
by measuring flow. 

The role of flow fluctuations and non-flow effects is another long standing problem 
that received a lot of attention and significant progress has been made in the recent 
couple years. In particular, the role of fluctuations in the initial system geometry [21] 
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Figure 4: Pion V2(pt) in two scenarios [19J. 
Solid line indicates the results using 
parameters best fit to spectra. 



defined by nuclear participants (25] (interacting nucleons or quarks) has been greatly 
clarified [2SJ E7J 12H1 EHJ EO]- The following picture emerges: Anisotropic flow is defined 
as the correlations to the reaction plane spanned by the impact parameter and the 
beam axis. At fixed impact parameter, the geometry of the participant zone fluctuates, 
both, in terms of the value of the eccentricity as well as the orientation of the major 
axes. Then the anisotropy develops along the plane spanned by the minor axis of the 
participant zone and the beam direction, the so called participant plane. As the true 
reaction plane is not known and the event plane is estimated from the particle azimuthal 
distribution "defined" by the participant plane, the apparent (participant plane) flow 
appears to be always bigger (and always "in-plane", V2,pp > 0) compared to the "true" 
flow as projected onto the reaction plane (see Fig. [7]). 

It was noticed in [25] that in collisions of heavy nuclei the fluctuations in the 
eccentricity (e x ,e y ) = {{{o~y — o%)/(o% + ^)), (faly/i&y + a x))) can be well described 
by two-dimensional Gaussian. What is not trivial is that for such Gaussian fluctuations 
the higher cumulant flow (v{n},n > 4) is not only insensitive to non-flow but also to 
eccentricity fluctuations. All of higher cumulants are exactly equal to the "true" flow, 
namely as given by projection onto the reaction plane. At the same time, the apparent 
(participant plane) flow become unmeasurable in a sense that flow fluctuations could 
not be separated from non-flow contributions by means of correlation measurements. 

An important conclusion from that study was that in most cases (except, probably, 
in mid-peripheral and peripheral Cu+Cu collisions, when the Gaussian approximation 
breaks [29J, and e part {4:} does not agree with e st d [26J ) the measurements of higher 
cumulant flow values provide the elliptic flow relative to the true reaction plane. 
Similarly, the same value is given by several other methods, such as Lee- Yang Zeroes, 
Bessel Transform, and fit to q-distributions. This greatly simplifies the comparison, 
e.g. of hydrodynamical calculations to the data, as it says that in such calculation one 
should not worry how to take into account the fluctuations in the initial eccentricity 
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(which is a non-trivial task) but just compare to the "right" measurement, e.g. t>2{4}. 
This understanding also allowed to appreciate some earlier calculations with uRQMD 
model [311 E2]. There, it was shown that using higher cumulants and/or LYZ method 
one indeed can measure the elliptic flow very well, but it was not at all clear why there 
have been observed no traces of the effects of flow fluctuations, which were expected in 
this model. 

Unfortunately this progress in understanding the nature of fluctuations does not 
help in resolving the problem of measuring flow fluctuations (in the participant plane) 
and non-flow. Strictly speaking, to make any estimates of those one is required to make 
assumptions. Most often to suppress non-flow contribution the azimuthal correlations 
between particles with large rapidity separation are used. The problem with this method 
is that there is no reliable estimates of how well it suppresses non flow and also how 
much the flow fluctuations (in this case correlations) change after imposing such a cut. 

At this conference the PHOBOS [33] and the STAR [M] collaborations presented 
their revised (compared to QM'06) results on flow fluctuations. These results are in good 
qualitative and quantitative agreement, see Figs. |5j[6j In [34J a conservative approach is 
taken and only upper limits on fluctuations are reported. The PHOBOS Collaboration 
uses estimates of non-flow effects from correlations with large rapidity separations and 
report more restrictive range for fluctuations. Both agree that the current measurements 
exhaust the (nucleon) eccentricity values obtained in MC Glauber model and in this 
sense somewhat favor models which predict smaller relative fluctuations, such as the 
CGC model or MC Glauber taking into account constituent quark substructure. 

Another important and interesting direction that is just started to be explored is 
the role of the non-zero initial flow velocity profile. The first obvious candidate here is 
to take into account the initial velocity gradient along the impact parameter, Fig. [SJ As 
shown in [35] such a gradient directly contributes to the in-plane expansion rate (see 
Eq. 22 in [35]). I would draw attention to the fact that those effects naturally also lead 
to directed flow (see the same Eq. 22), which too briefly addressed in |36j. It will be 
very interesting to compare the calculations in such a model to a very precise data from 
STAR [37] on directed flow obtained with the reaction plane determined by neutrons in 
the Zero Degree Calorimeter. The relation to other models [Ml EH] predicting non-trivial 
dependence of directed flow on rapidity would be also very interesting. Speculating on 
this subject one would also notice that viscous effects must also play an important role 
in such a scenario. 

It was shown in [40J that the effect of the strong parity violation that lead to non 
equal number of right and left fermions (quarks) in the presence of strong magnetic 
fields of colliding nuclei would result in charge separation (preferential emission of same 
charge particles) in the direction perpendicular to the reaction plane. Such anisotropy, 
which very much resembles "out-of-plane directed flow" can be addressed with the help 
of three-particle correlations [JT] by measuring (cos(0 a + (fip — 2^ rp)), where Q)/ g are 
azimuthal angles of two (same or opposite) charged particles, and ^pp is the reaction 
plane angle. The estimates [42J (see also talk by H. Warringa at this conference) 
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Figure 6: Elliptic flow fluctuations as 
estimated by PHOBOS Collaboration. 
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Figure 11: "Collapse" of directed flow as 
discussed in [35]. 



indicate that the effect is strong enough to be observed in heavy ion collisions. At this 
conference the STAR Collaboration reported the preliminary results [33], see Fig. [9j 
that qualitatively agree with estimates presented in [00], 42J. Note, that the correlator 
(cos(0 a + 0^ — 2^/ hp)) is P-even and contain contributions from other effects not related 
to parity violation. A careful analysis of such a contribution is obviously needed before 
any strong conclusion can be drawn from these measurements. 

Coming years promise many new interesting data from low energy RHIC run and, 
of course, from LHC. The main interest in the low energy RHIC scan, anisotropic flow is 
no exception, is the search for the QCD critical point. The scan would cover the energy 
region from top AGS energies, over the CERN SPS, and higher. In terms of anisotropic 
flow two major observables to watch would be a possible "wiggle" in v-ije dependence on 
particle density [33] and "collapse" of directed flow [32] ■ RHIC also has plans to extend 
its reach in terms of energy density using uranium beams. From the first estimates and 
ideas of using uranium beam we now have real detailed simulations [3S] of such collisions 
with developed methods for a selection of desired geometry of the collision. 

The predictions for the LHC are rather uncertain, though most agree that the 
elliptic flow will continue to increase [47J, partially due relatively smaller contribution 
of viscous effects. Simple extrapolations [48, 49j of the collision energy dependence of 
t>2 look rather "reliable'. Note that there exist calculations predicting decrease of the 
elliptic flow [2)]. Another important observation is an increase in mass dependence 
(splitting) of v 2 (pt) due to a strong increase of radial flow. 

In summary, we have had very exciting years of anisotropic flow study, which greatly 
enriched our understanding of ultra-relativistic nuclear collisions and multiparticle 
production in general. We are looking forward for new physics with LHC and RHIC. 

I thank A. Poskanzer, R. Snellings, and A. Tang for numerous fruitful discussions. 
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